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VISUAL  G INFRARED  ATMOSPHERIC  TRANSMITTANCE 
VISUAL  6 INFRARED  ATMOSPHERIC  SCATTERING  AND  ABSORPTION 
SURFACE  METEOROLOGY  AND  ATMOSPHERIC  RADIOMETRIC  TRANSFER 
IN  SITU  ATMOSPHERIC  MEASUREMENTS 

GROUND-BASED  RADIOMETRY 

Thii  a.  v 

Eastman  Kedak  Compaimrecently  completed  % contract  with  the  United  States 
Air  Force  Air  Weather  Service  to  conduct  research  into  the  transfer  of 
atmospheric  energy  in  the  visible  and  near-infrared  portions  of  the  elec- 
tromagnetic spectrum.  Its  objective  was  to  produce  a spectral  model  of 
path  transmittance,  radiance,  and  ground-level  irradiance  that  could  be  re 
lated  to  meteorological  observations  through  the  use  of  simultaneous,  in- 
situ  data  collections.  — - 


AStiTlQATlOjjl 


hcuwitv  ct*l»i»iC*Tiow  or  Twit  Moiwm  Of  hn«< 


The  instrumentation  package  developed  under  subcontract  to  Eastman  Kodak 
by  the  Albuquerque  Division  of  EGftG  (aee~Vu!  1)  was  used  in  a 14-month  data- 
collection  program.  The  program  had  two  objectives:  (1)  To  collect  spectral 
radiometric  data  on  the  atmosphere  and  on  ground- level  surface  irradiance 
under  the  broadest  possible  range  of  climatic/meteorological  conditions,  and 
(2)  to  collect  radiometric  and  meteorological  data  in  such  a way  that  would 
permit  them  to  be  related  mathematically.  In  this  field  collection  portion 
of  the  task,  data  was  obtained  for  18  different  climatic  conditions.  The 
collection  program,  which  lasted  18  months  and  covered  22,0^0  miles  of  the 
continental  United  Stqies,  required  extensive  coordination  with  other  organi- 
zations to  establish  acceptable  test  sites,  accumulate  correlated  data,  and 
conduct  concurrent  experiments  with  other  researchers  in  the  field  of  atmo- 
spheric sciences. 


^This  volume  (2)  of  the  final  report  describes  the  design  of  the  data  collec- 
tion program  conducted  in  the  field,  the  ranges  of  important  experimental 
parameters,  the  measurements  and  supportive  observations  that  were  made,  and 
the  analytical  software  developed  to  perform  the  statistical  analyses  on  the 
resulting  data. 
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SECTION  I 


INTRODUCTION 


The  instrumentation  package  developed  under  subcontract  to  Eastman 
Kodak  by  the  Albuquerque  Division  of  EGfjG  (see  Vol  1)  was  used  in 
a 14-month  data-col lection  program.  The  program  had  two  objectives: 
(1)  To  collect  spectral  radiometric  data  on  the  atmosphere  and  on 
ground- level  surface  irradiance  under  the  broadest  possible  range  of 
climatic/meteorological  conditions,  and  (2)  to  collect  radiometric 
and  meteorological  data  in  such  a way  that  would  permit  them  to  be 
related  mathematically. 


This  volume  (2)  of  the  final  report  describes  the  design  of  the  data 
collection  program  conducted  in  the  field,  the  ranges  of  important 
experimental  parameters,  the  measurements  and  supportive  observations 
that  were  made,  and  the  analytical  software  developed  to  perform  the 
statistical  analyses  on  the  resulting  data. 


SECTION  2 


ITINERARY  FOR  THE  COLLECTION  PROGRAM 

2.1  ROUTE  PLANNING  BASED  ON  REGIONAL.  CLIMATOLOGY 

To  make  spectroradiometric  measurements  under  all  the  meteorological 
conditions  that  were  to  be  mathematically  modeled  required  careful  plan- 
ing of  the  itinerary.  Data  collection  was  limited  geographically  to  the 
continental  U.S.,  and  measurements  had  to  be  completed  within  approximately 
one  year.  Fortunately,  quite  diverse  climatic  conditions  could  be  sampled 
in  the  U.S.by  means  of  proper  timing  with  respect  to  seasons.  Several 
climatological  descriptions  of  the  U.S.  were  found  that  were  utilized  in 
developing  the  collection  itinerary.  These  descriptions  indicated  that 
the  climate  over  large  contiguous  areas  is  somewhat  similar,  and  that  the 
climate  in  even  widely  separated  areas  sometimes  shows  like  characteristics. 
Because  nature  displays  no  sharp  distinctions,  overlapping  characteristics 
and  transition  zones  were  frequently  observed.  Although  the  positions  of 
the  transition  zones  are  not  fixed  permanently,  they  tend  to  be  relatively 
stationary.  . 

*,  5 

The  two  basic  methods  for  classifying  climate,  the  empiric  and  the  genetic, 
were  used  on  this  project.  Empiric  classifications  are  based  on  (1)  observ- 
able climatic  elements  such  as  temperature,  humidity,  and  precipitation,  and 
on  (2)  the  effects  of  these  elements  in  the  form  of  vegetation  zonation. 
Genetic  classification,  on  the  other  hand,  is  based  on  dynamic  weather  forces 
such  as  air  mass  type,  and  its  origin  and  the  nature  of  its  motion.  Although 
the  empiric  classification  was  used  to  develop  the  itinerary,  the  genetic 
classification  was  used  to  describe  the  air  system  sampled. 

2.2  THE  U.S.  ARMY  CLIMATIC  CLASSIFICATION  OF  THE  U.S. 

The  two  most  commonly  used  elements  in  empiric  classifications  are  precipita- 
tion and  temperature.  However,  temperature  alone  is  unsatisfactory  as  a basis 
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for  classification  because  it  would  provide  no  distinction  between  humid  and 
desert  regions.  Mean  precipitation,  although  an  equally  natural  element 
for  climatic  tone  classification,  is  not  necessarily  an  accurate  measure 
because  of  the  implications  of  seasonal  distribution  and  evaporation.  For 
example,  the  amount  of  precipitation  that  produces  high  humidity  conditions 
in  cold  temperatures  may  well  produce  dry  deserts  in  tropical  climates. 

As  a result,  a climate  classification  system  to  be  useful,  must  consider 
both  temperature  and  precipitation.  According  to  Visher,1  the  U.S.  Army 
developed  such  a method  in  1943  that  distinguished  twelve  zones  with  respect 
to  these  two  elements.  Figures  2-1A  through  2-11  show  the  locations  of  the 
zones  for  each  month  of  the  year.  For  this  project,  data  was  obtained  from 
as  many  of  the  12  zones  as  possible.  Although  many  of  the  zones  exist  at  the 
same  time,  the  areas  they  encompass  vary,  and  as  stated  above,  the  transition 
boundary  cannot  be  accurately  predicted.  For  this  reason,  it  was  planned  to 
acquire  data  in  a particular  zone  when  that  zone  was  large  and  from  a site 
that  did  not  lie  too  near  a transition  boundary.  In  addition,  the  travel 
route  was  laid  out  with  the  idea  of  minimizing  the  travel  time  between  acqui- 
sition sites  of  the  zones. 

The  number  of  days  allocated  to  each  site  was  based  on  the  probability  of 
obtaining  acceptable  weather  during  the  anticipated  period  of  visit  as  indi- 
cated in  publications  giving  climatological  data.2  Weather  was  classified  as 
acceptable  when  sky  conditions  were  from  clear  to  0.3  cloud  cover.  Rather  than 
"chase  the  weather,"  it  was  decided  in  advance  to  limit  the  stay  at  each  site 
to  the  number  of  days  that  had  been  allocated  whether  the  particular  meteor- 
ological/climatic condition  had  prevailed  or  not. 

By  using  the  monthly  climatology  maps  and  estimating  the  number  of  days  re- 
quired at  each  site,  a map  of  the  desirable  locations  was  prepared.  The  use 
of  similar  maps  of  military  bases  and  Rawindsonde  launch  sites  as  overlays  for 


1 S.  Visher,  Climatic  Atlas  of  the  United  States,  Harvard  University 
Press,  Cambridge,  Mass,  1954,  pp  373-376. 

2 H.  Conway  and  L.  Liston  (Editors),  The  Weather  Handbook,  Conway  Research 
Inc,  Atlanta,  Ga,  1974. 
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Q malic  lofKt  m Jenuery 
(Army  cleivihcetion  of 
l*4J):  1,  very  hot,  dry  (tem- 
perature over  46*F,  precipi- 
tation lea  then  4 in.) ; J , hoc, 
dry  (d»*-46\  lea  then  I 
in.);  J,  hut,  hunud  (4t*-46*, 
over  ) in  );  4.  warm.  dry 
(lO'-d**,  leu  then  2 in  );  J, 
Verm,  humid  (50*-6*\  2-4 
in.);  6,  mild,  dry  ()2*-S0*. 
lea  then  I in.);  7, mild, 
humid  ( >2*— JO*.  1-5  in.);  I, 
mild,  wet  (12  *-50*.  over  5 
in  );  f.  cold,  dry  (I4*-J7*. 
lea  then  I in.);  10.  cold,  hu- 
mid (>«*-»*.  l-J  in.);  II, 
cold,  wet  (l4*->2\  over 
J in.);  12.  very  cold  (be- 
low 14*). 


Qimeric  tones 
in  Februery 
(Army  deaificKion). 


Oimetic  nnee 
in  March 

(Army  daaification). 


Figure  2-1.  U.S.  Army  Classification  of  United  States 
Climate  by  Months 
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this  map  made  it  possible  to  determine  the  final  group  of  potential  sites. 

These  sites  are  listed  by  month  in  Table  2-1  in  the  order  dictated  by  the 
climological  study. 

Due  to  the  requirement  to  obtain  upper  air  soundings  at  each  acquisition 
location,  it  was  necessary  to  set  the  equipment  up  within  25  n.  mi.  of  the 
national  RAOB  launch  sites.3  Figure  2-2  indicates  the  locations  of  these 
facilities  by  attached  labels.  To  satisfy  power  requirements  and  physical 
security,  the  primary  sites  were  located  at  military  installations1'  wherever 
possible.  Many  of  these  turned  out  to  be  coincident  with  the  RAOB  launch 
sites.  Because  go  or  no-go  conditions  for  data  acquisition  and  hence  RAOB 
launches  were  often  not  known  until  shortly  before  launch  time,  the  setups 
adjacent  to  RAOB  sites  facilitated  decision  making  and  communication  with 
weather  personnel. 

2.3  BASE  LOGISTICS 

Overall  coordination  for  setup  of  sites  at  the  bases  was  provided  by  the  AIR 

. > 

Systems  Command/LGYJ . Each  potential  data  collection  site  required  advance 
< arrangements  with  respect  to  base  entry,  a suitable  site  for  trailer-setup, 

power,  sleeping  and  eating  facilities  for  a two-man  crew,  and  physical  security. 
These  arrangements  were  made  by  LGYJ  personnel  by  telephone,  and  Air  Force  con- 
tractor travel  orderfc  were  provided.  Throughout  the  collection  phase,  the 
code  name  SENIOR  VICTOR  was  used  to  identify  this  project.  Table  2-2  lists  the 
locations  of  sites  used  during  this  project. 

2.4  AIR  WEATHER  LOGISTICS 

N Of  major  importance  to  the  project  was  the  ability  to  obtain  accurate  and 

frequent  meteorological  observations  at  or  near  the' base  of  operations.  These 

3 NDAA  Weather  Service  Offices  and  Stations,  Jan  ^75,  and  U.S.  Dept  of 
Commerce  RAOB -RAW IN  Altitude  Summary,  Dec  1974. 

• Air  Force  Magazine,  "Air  Force  Magazine's  Guide  to  U.S.A.F.  Bases  at 

< Home  and  Abroad",  May  1975,  pp  142-154. 
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LIST  OF  SITES  BY  MONTH  PLANNED  FOR  THE  COLLECTION  PROGRAM 
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23A  Dulles  A/P,  Va.  400  Vint  Hill  A ray  Depot  7 Amy  3B47N/7741N 

24  Dayton,  Ohio  400  Wright-Patterson  AFB  12  AFLC.AFSC  29SON/8403K 

25  Qaaha,  Neb.  650  Offut  AFB  7 SAC  4107N/9SSS* 


LIST  OF  SITES  BY  MONTH  PLANNED  FOR  THE  COLLECTION  PROGRAM 


ometric  Laboratory  Route  and  Anticipated  Measurement 


TABLE  2-2 

ACTUAL  SITES  USED  IN  THE  DATA  COLLECTION  PROGRAM 
(Project  SENIOR  VICTOR) 


were  of  three  types:  normal  synoptic  surface  weather  observations,  Rawindsonde 
observations,  and  reconstructed  imagery  from  meteological  satellites. 

2.4.1  Synoptic  Data. 

Although  the  trailer  was  equipped  with  a meteorological  system,  data  was  ob- 
tained from  the  weather  service  at  each  site  of  operation.  When  there  was 
serious  disagreement  between  the  trailer  system  and  the  weather  service,  the 
latter  was  the  one  used  in  data  analyses.  The  surface  observations  made  by  the 
trailer  weather  system  for  each  site-day  are  summarized  in  Appendix  A of  this 
report . 

The  surface  data  collected  by  the  on-board  weather  system  was  summarized  by  a 
computer  program,  SUNSHE,  that  plotted  and  computed  various  parameters  as  a 
function  of  Greenwich  Mean  Time  (GMT).  To  provide  an  unambiguous  time  base, 

GMT  was  used  throughout  the  collection  program.  Relative  humidity  (RH)  and  tem- 
perature (T)  were  reduced  to  surface  vapor  pressure  (VP)  by  substituting  them 
in  the  following  equation: 

VP(t)  - RH(t)  x [1.186  ♦ 4.514  x 10"2T(t) 

♦ 3.921  x 10'3T2(t)  - 3.456  x 10'5T3(t) 

♦ 5.548  x 10'7T4(t)J 

Rawindsonde  data  was  available  three  times  each  day:  at  1200Z,  1800Z,  and 
2400Z;  and  vapor  pressure  was  interpolated  for  these  times.  The  pressure  and 
wind  data  were  reported  over  the  entire  day  of  collection. 

2.4.2  Rawindsonde  Data 

By  prior  arrangement  with  the  weather  service,  one  additional  Rawindsonde  balloon 
was  launched  at  mid-day  on  days  of  trailer  operation.  In  general,  the  on-site 
coordination  for  these  was  performed  by  the  field  crews  with  advance  notification 
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of  the  launch  requirement  arranged  through  AFAWS.  A list  of  the  Rawindsonde 
Atmospheric  Observation  Balloon  (RAOB)  stations  used  during  this  project  can 
be  found  in  Table  2-2. 

RAOB  station  data  was  transmitted  to  AFAWS  for  computer  processing.  The  data 
when  finally  received  listed  precipitable  moisture,  wind  data,  temperature, 
humidity,  density,  and  pressure,  as  a function  of  altitude  above  the  local 
station  elevation.  An  example  of  these  data  is  shown  in  Figure  2-3  for  Denver, 
Colorado,  21  October  1975  at  18002.  The  surface  pressure  and  integrated  precipitable 
moisture  scale  height,  which  were  the  key  parameters  in  atmospheric  radiometric 
data  scaling,  are  explained  in  detail  in  Section  2 of  Volume  3 of  this  report. 

2.4.3  Meteorological  Satellite  (Metsat)  Imagery 

AWS  also  provided  paper  print  copies  of  the  reconstructed  imagery  from  the  digital , 
meteorological  satellite  signals  in  the  visible  spectral  band.  An  example  is  shown 
in  Figure  2-4.  The  imagery  was  used  periodically  to  diagnose  cloud  patterns  or 
atmospheric  conditions  for  the  macro-scale  near  the  site  of  operation.  Although 
sometimes  useful  in  understanding  large  patterns  of  weather, the  scale  was  too 
small  to  be  of  practical  value  in  the  project. 


Figure  2-3.  Exaaple  of  RAOB  Data  for  Denver,  Colorado, 
on  21  October  1975,  1800Z 
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SECTION  3 


SUMMARY  OF  DATA  COLLECTIONS 

3.1  OVERVIEW  OF  THE  EXPERIMENT 

The  data  collection  program  was  conducted  from  August  1975  through  October 
1976,  and  the  mobile  laboratory  visited  62  sites  within  the  United  States. 
Data  collections  were  made  both  before  and  after  local  true  noon  so  that 
each  day-site  could  potentially  yield  two  separate  collections.  As  this  oc- 
curred 27  times  during  the  program,  the  total  number  of  data  collections  was 
89  (2x27+35).  Appendix  A of  this  volume  summarizes  these  collections,  list- 
ing the  exact  values  or  ranges  of  physical  parameters  that  were  important  to 
this  experiment.  This  summary  is  also  shown  graphically  in  Figure  3-1. 

The  column  CONF  COEF  (confidence  coefficients)  in  Figure  3-1  was  developed 
to  separate  data  collections  into  (1)  those  that  would  be  usable  in  the 
analysis  phase  of  the  project,  (2)  those  that  would  be  usable  as  test  cases, 
and  (3)  those  that  would  not  be  usable  at  all.  Since  all  radiometric  data 
measurements  were  made  concurrently  with  spectral  transmittance  measurements, 
the  coefficient  was  computed  on  the  basis  of  the  range  of  air  mass  and  the 
number  of  samples  corresponding  to  the  transmissometer  measurements.  Spe- 
cifically, the  coefficient  was  the  product  of  the  air-mass  range  and  the 
square  root  of  the  sample  size,  as  shown  below,  and  reflected  the  equations 
associated  with  the  standard  error  of  regression  as  discussed  by  Rickmers5 
and  Todd. 

Confidence  Coefficient  = AM*  x t/""N 

Where,  AM  is  the  range  of  the  air  mass,  and  N is  the  number  of  spectral 
samples . 


5 A.  D.  Richmers  and  H.  N.  Todd,  Statistics,  an  Introduction,  McGraw  Hill 
Inc,  New  York,  N.Y,  1967,  p.258. 

* AM  is  an  approximation  of  the  quantity  tQ  °y.x  discussed  in  reference  5. 

In  addition,  it  is  assumed  that  the  prediction  of  y is  for  the  average 
value  of  x. 
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Four  ranges  of  the  confidence  coefficient  were  used  to  separate  data  into 
excellent,  good,  fair,  and  poor  categories.  Of  the  89  total  collections, 

30  were  considered  excellent,  9 good,  13  fair,  and  37  poor.  Only  excellent 
and  good  collections  were  used  in  formulating  the  final  atmospheric  model; 
the  fair  and  poor  data  were  only  used  as  test  cases  or  not  at  all.  In 
general,  the  fair  and  poor  ratings  represented  an  insufficient  range  or 
sample  size.  They  typically  resulted  from  non-homogeneous  cloud  cover  that 
prematurely  terminated  measurements;  however,  a few  were  the  result  of 
instrument  failures. 

3.2  DISTRIBUTION  OF  SAMPLES  BY  PARAMETER 

Figures  3-2  through  3-5  are  histograms  of  the  number  of  data  collection 
spectral  samples  for  each  physical  parameter  and  for  the  entire  experiment. 

Each  point  of  the  histograms  represents  a set  of  four  simultaneous  spectral 
measurements  from  each  of  four  spectroradiometers . Only  data  sets  having 
confidence  coefficients  in  the  good  or  excellent  categories  are  included  in 
the  histograms.  Although  a good  range  of  the  parameters  is  evidenced,  the 
samples  did  not  represent  a uniform  collection  distribution  (e.g.,  solar 
altitude  is  strongly  peaked  at  30°  and  skewed  toward  the  high  solar  altitudes) . 
In  formulating  the  equations  of  the  atmospheric  model,  this  artifact  of  data 
collection  was  compensated  for  as  follows.  The  equations  resulting  from 
regression  analysis  were  adjusted  on  the  basis  of  examination  of  systematic 
sampling  of  data  sets.  Table  3-1  summarizes  the  important  statistics  of  the 
histograms . 

TABLE  3-1 

SUMMARY  OF  STATISTICS  CORRESPONDING  TO  HISTOGRAM 
FIGURES  3-2  THROUGH  3-5 


Parameter 

Mean 

Lowest 

Highest 

Solar  Altitude  (degrees) 

36.3° 

O 

00 

iH 

1 

80.8° 

Moisture  Scale  Height  (cm) 

1.61  cm 

0.22  cm 

4.58  c 

Terrain  Reflectance  (§  650  nm) 

17.6% 

4.0% 

60.3% 

Terrain  Altitude  (feet) 

2127  ft 

19  ft 

5400  ft 

3-5 


Mean  * 36.3 

H2  ‘ Std.  Dev.  = 19.3 

N = 947 


Figure  3-2.  Histogram  of  Number  of  Spectral  Samples  by  Solar  Altitude 


Std.  Dev 


Figure  3-3.  Histogram  of  Number  of  Spectral  Samples  by  Moisture  Scale  Height 


Figure  3-4.  Histogram  of  Number  of  Spectral  Samples  by  Terrain  Reflectance  (650  nm) 


Histogram  of  Number  of  Spectral  Samples  by  Terrain  Altitude 


SECTION  4 

MEASUREMENT  TYPES  AND  SEQUENCES 


4.1  DESCRIPTION  OF  MEASUREMENTS 

The  design  of  the  collection  program  included  a description  of  the  types  fincl 
frequency  of  certain  measurements  on  each  site.  This  plan  was  implemented  by 
the  field  operators  and  used  to  identify  the  measurements  of  a common  type 
in  later  analyses.  A number  of  priorities  was  also  established  as  a guide 
for  making  operational  decisions  in  the  field,  as  follows: 

Priority  1:  Measurements  of  greatest  significance  that  were 
to  be  made  routinely  each  day  of  operation. 

Priority  2:  Measurements  required  at  each  site,  but  not  neces- 
sarily each  day.  These  were  to  be  made  at  time 
intervals  not  interfering  with  Priority  1 collec- 
tions and  on  days  following  successful  Priority  1 
collections. 

Priority  3:  Measurements  of  a specialized  nature  that  were  to 
be  accomplished  as  specially  requested  operations. 

They  were  not  to  be  made  unless  Priority  1 and  2 
collections  had  already  been  satisfied  for  each  site. 

The  codes  assigned  for  measurements  and  a brief  description  of  the  wavelength 
and  orientation  limitations  corresponding  to  each  is  contained  in  Table  4-1. 

It  was  necessary  to  define  the  exact  sequence  of  measurements  so  that  collec- 
tion of  the  most  critical  measurements  could  be  made  and  still  provide  time  for 
alternate  observations.  For  this  purpose,  10  unique  instrument  configurations 
were  defined.  Table  4-2  lists  the  operation  sequences  in  terms  of  the  servo 
positions  and  special  geometry  requirements.  Since  there  were  two  identical 
irradiometers,  selection  of  the  instrument  to  make  the  irradiance  measurement 
was  left  to  the  operator's  discretion.  A predetermined  sequence  of  92  angle 
pairs  were  used  for  the  azimuth  and  elevation  positions  of  the  sky  radiometer 
when  the  "sky  map"  measurements  were  made.  The  angles  were  selected  to  provide 
complete  sky  hemisphere  coverage  with  minimum  overlap  of  the  1S°  field-of-view 
optics.  Table  4-3  lists  the  angular  positions  in  the  order  of  operation.  All 
angular  measurements  are  relative  to  the  sun's  azimuth  (clockwise  rotation, 
from  top)  and  the  true  instrument  horizontal. 
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DESCRIPTION  OF  RADIOMETRIC  MEASUREMENT  PROGRAM 


Rawindsonde  Upper  Air  1 Local  RAOB  launch  task  at  or  near  local  noon.  No  field 

Sounding  data  required. 


TABLE  4-2 


OPERATING  SEQUENCES  (CONT'D) 


aqoisition 

NO. 

IRRADIOKETER 

IRE 

IRRADIOMETER 

IR 

SKY  RADIOMETER 

TRANSMI SS0YETE0 

3 

Atlauth 

Elevation 

Dice  Block 

Wavelength 

Description 

4*0 

Aso-180* 

Aso-180* 

4*0 

90* 

0* 

90*-2SE  or  90* -SE 

ELo 

Tm 

— 

mmm 

— 

350-1200,10 

350-1200,10 

350-1200,10 

350-1200,10 

Boris.  Mil  Sky 

Point  Mode  Sky 

Solar  Disc 

Xrradlance 

Irradlanee 

Radi  a nee 

Radiance 

6 

Atlauth 

novation 

Disc  Block 

Wavelength 

Description 

4*0 

Aso-180* 

As,  Atc 

4*0 

90* 

0* 

D.S  Db 

do 

Tm 

MM 

Tes 

— 

1*50,650,850 

1*50,650,850 

1*50,650,850 

1*50,650,850 

Boris.  Mil  Sky 

Vert.  Back-Lit 

Mapping  Mode 

Solar  Dlse 

Irradlanee 

Radlanee 

7 

Atlauth 

Elevation 

Disc  Block 

Wavelength 

Description 

Aso-180* 

Asn-180* 

Asq-180* 

• 

4*0 

90* 

0* 

90*-2SE  or  90* -SE 

Ho 

• 

r T M 

350-1200,10 

350-1200,10 

350-1200,10 

350-1200,10 

Boris.  Daylleht 

Vert.  Back-Lit 

Point  Mode  Sky 

Solar  Dlse 

Irradlanee 

Irradlanee 

Radlanee 

Radlanee 

8 

Atlauth 

Elevation 

Dlae  Block 

Wavelength 

Description 

Aso-1800 

Aso-180* 

Asf  Ask 

4*0 

90* 

0* 

P,  P.F 

Eo 



aJU 

Ymm 

— 

1*50,650,850 

1*50,650,850 

1*50,650,650 

1*50,650,850 

■sum 

Solar  Dlse 

IffggHl 

Irradlanee 

Radlanee 

5 


TABLE  4-2 
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OPERATING  SEQUENCES  (CONT'D) 


AQDISITJOII 

HO. 


Aklauth 

Cwatlen 


WRAMOKETtR 
IRB 


XRRAD^HCTER 


SKY  RADIOMETER ] 

A •0-190* 


TRARSW  SSC'^T? 


Description 


TABLE  4-3 


SKY 

MAP 

MEASUREMENT 

ELEVATIONS 

AND  AZIMUTHS  IN 

THE  ORDER  OF 

OPERATION 

Point 

Point 

Point 

No. 

AZ 

EL 

No. 

AZ 

EL 

No. 

AZ 

I! 

2 

7 

7 

33 

122 

22 

64 

240 

67 

3 

7 

22 

34 

112 

7 

65 

240 

82 

4 

4 

37 

35 

127 

7 

66 

276 

67 

5 

22 

67 

36 

138 

22 

67 

279 

52 

6 

22 

52 

37 

138 

37 

68 

269 

37 

7 

22 

37 

38 

151 

52 

69 

252 

22 

8 

22 

22 

39 

154 

37 

70 

247 

7 

9 

22 

7 

40 

154 

22 

71 

262 

7 

10 

37 

7 

41 

142 

7 

72 

269 

22 

11 

40 

22 

42 

157 

7 

73 

277 

7 

12 

40 

37 

43 

172 

7 

74 

285 

22 

13 

48 

52 

44 

170 

22 

75 

288 

37 

14 

60 

67 

45 

174 

37 

76 

302 

22 

15 

73 

52 

46 

174 

52 

77 

292 

7 

16 

61 

37 

47 

168 

67 

78 

307 

7 

17 

56 

22 

48 

195 

37 

79 

318 

22 

18 

52 

7 

49 

187 

22 

80 

307 

37 

19 

67 

7 

50 

187 

7 

81 

307 

52 

20 

72 

22 

51 

202 

7 

82 

312 

67 

21 

80 

37 

52 

202 

22 

83 

330 

52 

22 

89 

22 

53 

212 

37 

84 

326 

37 

23 

82 

7 

54 

204 

52 

85 

334 

22 

24 

97 

7 

55 

204 

67 

86 

322 

7 

25 

105 

22 

56 

228 

52 

87 

337 

7 

26 

99 

37 

57 

231 

37 

88 

351 

7 

27 

99 

52 

58 

220 

22 

89 

351 

22 

28 

96 

67 

59 

217 

7 

90 

345 

37 

29 

120 

82 

60 

232 

7 

91 

356 

52 

30 

132 

67 

61 

236 

28 

92 

348 

67 

31 

125 

52 

62 

250 

37 

93 

360 

82 

32 

118 

37 

63 

253 

52 

4.2  DESCRIPTION  OF  DAILY  OPERATION 

The  four  MRMS  radiometric  instruments  (two  irradiometers,  a sky  radiometer, 
and  a solar  disc  radiometer)  permitted  as  many  as  four  measurements  to  be 
made  simultaneously,  such  measurements  being  considered  a measurement  set. 

Each  set  included  measurement  of  the  solar  disc  radiance,  which  is  used  in 
the  determination  of  atmospheric  transmittance. 

Each  of  the  four  MRMS  instruments  is  operated  by  command  software,  which 
determines  the  wavelengths  at  which  measurements  are  taken.  In  general,  all 
measurements  except  sky  radiance  in  the  map  mode  were  made  from  350  to  1200 
nanometers  in  10  nanometer  increments.  Because  of  time  constraints,  sky 
radiance  in  the  mapping  mode  was  only  measured  at  three  wavelengths:  450, 

650,  and  850  nanometers.  Other  measurements  of  the  same  set  had  to  be  made 
at  these  wavelengths  also. 

The  normal  sequence  of  operations  consisted  of  Priority  1 collections  inter- 
spersed, where  possible,  with  Priority  2 collections  according  to  the  schedule 
of  operations  1 through  10.  Sequence  1 through  10  of  Table  4-2  was  repeated 
several  times  as  daylight  hours  permitted.  Because  solar  altitude  changes 
occur  more  rapidly  near  dawn  and  dusk,  it  was  also  necessary  to  group  the 
instrument  sequences  into  solar  altitude  ranges.  For  this  purpose,  the  field 
crews  used  the  following  acquisition  sequences: 

Above  10°  Solar  Altitude:  Make  acquisitions  1 through  10. 

Between  5°  and  10°  Solar  Altitude:  Make  acquisitions  1 through  4; 
replicate  in  reverse  order  as  time  permits. 

Below  5°  Solar  Altitude 

On  Odd-Numbered  Days  of  Collection: 

Make  acquisition  9.  Repeat  first  wavelength  measurement 
after  spectral  scan  is  completed  to  serve  as  a monitor 
of  level  change  during  time  of  collection. 

Make  acquisition  10.  Use  horizontal  full  sky  irradiance 
measurement  to  serve  as  a monitor  of  sky  radiance  level 
change  during  time  of  collection. 
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On  Even  Numbered  Days  of  Collection: 


Make  acquisition  3.  Repeat  first  wavelength  after  spectral 
scan  is  completed  to  serve  as  a monitor  of  level  change 
during  time  of  collection. 

Make  acquisition  4.  Use  horizontal  full  sky  irradiance 
measurement  as  a monitor  of  sky  radiance  level  change 
during  time  of  collection. 

4.3  OPERATOR  COMMENT  INFORMATION 

The  measurements  discussed  above  were  continued  throughout  the  day  as  weather 
permitted,  and  the  data  was  recorded  automatically  on  tape  cassettes.  At 
the  end  of  an  acquisition  day,  supporting  data  for  each  acquisition  was  placed 
on  the  data  tapes  as  header  block  information.  The  supporting  data  for  each 
acquisition  was  entered  in  coded  form  (see  Figure  4-1)  and  included  the  type 
of  measurement  made  by  each  irradiometer,  cloud  information,  and  an  estimate 
of  data  reliability. 

Cloud-cover  data  consisted  of  location(s)  relative  to  the  sun,  degree  of 
cloud  cover,  and  the  cloud  type.  This  data,  based  on  observations  at  the 
time  of  the  data  acquisition,  was  designed  to  supplement  the  sky  camera  acqui- 
sitions, which  were  made  at  approximately  15-minute  intervals.  For  example, 
some  cirrus  clouds  that  could  go  undetected  in  the  sky  camera  acquisition 
might  significantly  affect  the  data. 

The  data  reliability  estimate  was  made  in  the  form  of  a "0"  to  "3"  rating  scale 
for  the  measurement.  A rating  of  3 meant  that  the  data  were  apparently  good, 
with  no  significant  changes  in  ambient  conditions,  no  instrument  problems, 
etc.  A "0"  rating  meant  that  data  should  be  rejected.  Data  rated  "1"  and 
"2"  had  certain  potential  problems  associated  with  them  that  required  judgment 
with  respect  to  their  usability. 

On-board  weather  measurements  were  acquired  with  every  radiometric  data  set. 
These  measurements  included  temperature,  relative  humidity,  barometric  pres- 
sure, surface  wind  speed,  and  direction.  Data  were  also  provided  by  the  Air 
Weather  Service  (AWS)  or  the  National  Weather  Service  (NWS)  at  most  sites. 

This  data  was  discussed  in  more  detail  earlier  in  this  report  (Section  3). 


ESTIMATE  OF  DATa  RELIABILITY : 

■3»  - Apparently  Good  Data  (No  significant  change  in  ambient 
conditions,  no  instrument  problems,  ate.) 

"2*  . Probably  Good  Data  (Notlcoabla  change*  in  ambient 
conditions  requiring  Judgments  in  classification, 
possibly  bad  data  points  requiring  correction,  etc.) 

"1"  - Suspect  Data  (Widely  varying  ambient  conditions, 

possibly  obscurstion  of  the  sun  during  measurements 
and/or  instrueMnt  unreliability.  Judgments  will 
have  to  be  made  on  the  usefullness  of  thia  data 
prior  to  incorporation  into  data  base.) 

"0"  - Bad  Data/No  Data  - Reject 


MEASUREMENT  CODE:  . 

"DH"  - Horisontal  Daylight 

"DV"  - Vertical  Daylight 

"SH"  - Horisontal  Skylight 

"SVF"  - Front  Lit  Vertical  Skylight 

■SVB"  - Back  Lit  Vertical  Skylight 

"NO"  - Instrument  Not  Operational 


CLOUD  COVER; 

Sky  conditions  during  the  data  acquisition. 

CLOUD  LOCATION  BY  SECTOR i 

Position  of  clouds  relative  to  the  sun  / 
(see  diagram)  as  observed  from  ground  L 
level.  T 

DEGREE  OF  CLOUD  COVER i 


Amount  of  cloud  cover  in  a specified  sector 
in  tenths. 

CLOUD  TYPE; 

"C"  - Cumulus 
"S"  - Stratus 
"R"  - Cirrus 


Figure  4-1.  Header  Block  Comment  Data 
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FORMAT  FOR  MEASUREMENT  COSMENTS  ON  MAGNETIC  TAPES 


A complete  set  of  panoramic  photographs  was  taken  of  the  surrounds  of  each 
measurement  set.  These  photographs,  taken  of  the  horizon  at  30-degree  inter 


vals  in  all  directions,  show  the  surround  terrain,  vertical  obstructions, 
and  any  physical  factors  that  might  influence  the  acquired  data.  In  addition, 
full-sky-hemisphere  photography  was  obtained  using  a "fish-eye"  lens  adapted 
to  the  same  camera.  These  images  were  taken  continually  throughout  the  mea- 
surement period.  Figure  4-2  is  an  enlargement  of  one  such  photograph.  When- 
ever it  was  practical,  aerial  photographs  were  taken  of  the  site  and  were 
used  to  determine  the  effective  albedo  reflectance  (see  Section  5 of  this 
report) . 
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SECTION  5 

SPECIAL  CONCURRENT  EXPERIMENTS 


5.1  DETERMINATION  OF  SITE  SURROUND  REFLECTANCE 

5.1.1  High  Altitude  Photographic  Coverage 

Under  the  direction  of  A.F./LGYJ,  U-2  flight  tests  were  scheduled  concur- 
rently with  the  atmospheric  measurements  made  by  the  mobile  laboratory. 

The  purpose  of  the  flights  was  to  obtain  color  imagery  that  could  be  used 
to  estimate  the  general  terrain  spectral  reflectance  of  the  site  surround 
terrain;  this  was  later  used  in  formulating  the  atmospheric  model.  The 
flight  tests  associated  with  the  atmospheric  measurements  are  listed  in 
Table  5-1.  Because  of  the  occasional  unavailability  of  ground  control  tar- 
gets, or  CORN  targets*  and  aircraft,  only  three  flights  were  actually  flown 
concurrently  with  atmospheric  measurements.  However,  most  flights  were  con- 
ducted within  2 to  10  days  of  the  atmospheric  measurements,  and  it  was  be- 
lieved that  the  time  span  between  flights  and  measurements  was  so  short  that 
no  appreciable  change  in  terrain  reflectance  would  occur.  One  flight  was 
actually  flown  one  year  later  but  at  nearly  the  same  time  of  the  month.  Be- 
cause considerable  distances  were  involved,  it  was  not  possible  to  overfly 
all  measurement  sites.  For  these,  terrain  reflectance  estimates  were  made 
using  a geobotanical  map  of  the  area. 

For  these  tests,  U-2  aircraft  were  used  flying  at  an  altitude  of  60,000  feet 
with  a pair  of  24-inch  panoramic  cameras.  The  film  load  consisted  of  SO-242 
High  Definition  Color  Film  in  the  forward  camera  and  SO-217  High  Definition 
Black-and-White  Film  in  the  aft  camera.  The  aperture  was  set  at  f/4.6  for 
the  forward  camera  and  f/3.5  for  the  aft.  The  exposure  time  was  varied  over 
a 6-ms  to  8-ms  range  for  the  test  depending  on  the  solar  altitude.  The  time 
of  photography  was  scheduled  with  the  midpoint  of  the  flight  to  occur  at  true 
noon.  This  eliminated  any  substantial  changes  in  solar  altitude  that  would 
have  necessitated  a change  in  exposure  during  the  test. 


* control  Range  Network,  prepared  and  displayed  by  Mead  Technologies 
Laboratories,  Dayton,  Ohio. 
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TABLE  5-1 


SUMMARY  OF  U-2  OVERFLIGHTS 


Field  Correspnding  Atmospheric 

No.  Date  Flown  Location  Measurement  Date(s) 


76-1 

27  Feb  76 

St.  Cloud,  Minn 

21-27  Feb  76 

76-3 

1 Apr  76 

Quillayoute,  Wash. 

3-6  Apr  76 

76-4 

21  Apr  76 

Moffet  Field,  Cal. 

13-26  Apr  76 

76-6 

4 Jun  76 

San  Diego,  Cal. 

12  May-1  Jun  76 

76-7 

9 Jun  76 

Tucson,  Ariz. 

7-17  Jun  76 

76-8 

22  Jun  76 

Winslow,  Ariz. 

28  Jun  76 

76-9 

29  Sep  76 

Albuquerque,  N.M. 

12-19  Sep  76 

76-10 

5 Oct  76 

Ogden,  Utah 

27-29  Sep  76 

76-11 

13  Oct  76 

Denver,  Colo. 

21-22  Oct  75 

CORN  targets  were  deployed  at  each  site  for  calibration  and  analysis  of 
reflectance.  Each  deployment  consisted  of  a 5-step  gray  scale,  a 3-step 
gray  scale,  a 12-  (or  6-)  step  color  target,  and  a tri-bar  resolution  target. 
On  some  tests,  not  all  the  targets  were  displayed.  The  smallest  display  con- 
sisted of  a 6-step  color  target  and  a 5-step  gray  scale  target.  Figure  5-1 
is  an  enlargment  made  from  high  altitude  color  photography  of  a typical  CORN 
deployment  at  a measurement  site.  On  one  flight  (76-4),  the  targets  were  dis- 
played at  a distance  of  20  miles  from  the  mobile  laboratory  location;  however, 
the  terrain  was  quite  similar  and  it  was  still  possible  to  use  the  targets  for 
calibration.  Four  passes  over  the  target  area  were  made  by  the  U-2  aircraft 
to  permit  enough  repetitive  photography  to  be  acquired  to  provide  the  analysts 
with  reasonable  confidence  in  the  terrain  reflectances  derived. 

5.1.2  Analysis  of  Color  Imagery 

Microdensitometry  of  the  color  imagery  resulted  in  red,  green,  and  blue  den- 
sities that  corresponded  to  each  CORN  panel.  By  regressing  the  measured 
densities  from  each  photograph  against  the  known  ground  reflectances  of  each 
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panel  relative  to  each  layer  of  the  color  film,  it  was  possible  to  produce 
calibration  between  density  and  integrated  reflectance.  The  spectral  reflec- 
tance data  was  integrated  numerically  with  respect  to  the  camera/ film  system 
spectral  response  to  produce  three  effective  reflections.  Large  areas  of  the 
surround  photography  were  then  scanned  with  a microdensitometer  to  estimate 
the  images'  average  red,  green,  and  blue  densities  corresponding  to  the  general 
terrain.  These  densities  were  then  passed  through  the  calibrations  to  obtain 
integrated  terrain  reflectances.  The  calibration  equations  were: 


log  Rr 

■ 

a0  * alDr  * a2°r  x Dg  * a3°r  x Db 

log  Rg 

m 

K * * b2Dg  X Dr  . bjDg  X Db 

l°g  Rjj 

s 

c0  * clDb  * =2Db  x Dr  * c3Db  x Dg 

Khare;  Dr,g.b 

are 

Status  A red,  green,  and  blue  densities  of  the 

CORN  panels, 

R . are  integrated  red,  green,  and  blue  reflectances  of  the 
r ,g,D 

CORN  panels,  and 

an  b.  c.  are  regression  coeficients. 

Table  5-2  lists  the  resulting  integrated*  reflectances  corresponding  to  the 
flight  tests  in  Table  5-1. 


* Integrated  relative  to  40  degrees  solar  altitude  daylight  irradiance, 
24- inch  Petzval  lens,  and  SO-242  color  film. 
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TABLE  5-2 


INTEGRATED  TERRAIN  REFLECTANCES  ESTIMATED  FROM  U-2  COLOR  PHOTOGRAPHY 

Percent  Reflectance 


Site 

Flight  No. 

Red 

Green 

Blue 

St.  Cloud,  Minn 

76-1 

40.7 

40.7 

41.7 

Quillayoute,  Wash 

76-3 

3.8 

4.5 

4.1 

Alameda,  Calif 

76-4 

11.2 

11.8 

6.8 

San  Diego,  Calif 

76.6 

14.5 

12.3 

8.3 

Tucson,  Ariz 

76-7 

26.3 

20.4 

15.1 

Winslow,  Ariz 

76-8 

23.4 

15.1 

14.1 

Albuquerque,  NM 

76-9 

18.6 

13.8 

12.3 

Ogden,  Utah 

76-10 

7.6 

6.1 

4.3 

Denver,  Col 

76-11 

12.3 

10.2 

8.1 

Average 

14.4 

10.2 

8.1 

5.1.3  Reconstruction  of  Terrain  Spectral  Reflectances 

The  complete  spectral  reflectance  of  the  surrounding  terrain  was  reconstructed 
using  an  Eigenvector  technique.  Two  sets  of  spectral  reflectance  data,  one 
for  plants  and  trees,  and  the  other  for  sands  and  soils, were  analyzed  to  produce 
a set  of  characteristic  vectors  for  the  two  sample  spaces.  Each  spectral  re- 
flectance curve  of  the  sample  spaces  could  then  be  reconstructed  by  multiplying 
the  vectors  of  the  corresponding  space  by  the  coeficients  computed  for  that 
specific  sample  data,  as  follows: 

RS(X)  - Ms(X)  + Ax  • Vsl(A)  ♦ A2  • Vs2(X)  ♦ A3  • Vs3(X) 
or, 

Rp(X)  - Mp(A)  ♦ Bx  • Vpl(X)  ♦ B2  • Vp2(X)  ♦ Bj  . Vp3(X) 

Where:  R (X)  is  the  spectral  reflectance  of  any  specific  soil  or 
s sand  sample, 

R (X)  is  the  spectral  reflectance  of  any  specific  plant  or 
P tree  sample, 
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are  the  mean  row  vectors, 
are  the  chracteristic  vectors, 

are  coefficients  corresponding  to  specific 
samples  in  the  original  sample  sets. 

The  spectral  reflectance  data  was  also  digitally  integrated  relative  to  the 
three  bands  of  response  corresponding  to  the  color  film,  camera  lens,  and 
daylight  irradiance  used  in  the  terrain  reflectance  flight  tests. 

Using  linear  regression  techniques,  a transformation  was  found  between  the 
integrated  effective  reflectances  and  the  vector  coefficients;  that  is. 


Ms,P<» 
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Given  any  set  of  integrated  red,  green,  and  blue  reflectances  (R^R^R^), 
it  was  possible  to  estimate  a spectral  reflectance  curve  for  soils  and 
plants  that  could  have  produced  the  measured  reflectances.  If  the  operating 
site  terrain  had  been  all  sand  or  all  forest,  then  the  estimated  spectral  re- 
flectance would  only  have  been  that  reconstructed  for  the  corresponding  spec- 
tral- data  space.  Because  this  was  not  generally  the  case,  a mixture  of  ground 
types  was  present.  From  visual  inspection  of  the  U-2  photography,  it  was 
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possible  to  estimate  the  fractions  corresponding  to  sands/soils,  and  to 
plants/trees.  In  addition,  an  assessment  of  the  fraction  of  water  and/or 
snow  surround  was  also  made.  These  four  surround  types  were  then  combined 
by  weighting  the  spectral  reflectance  corresponding  to  each  by  the  visual 
fraction  to  produce  an  overall  estimate  of  the  local  earth  albedo  reflec- 
tance. Water  and  snow  were  handled  as  constant  spectral  reflectances  and 
not  reconstructed  from  any  vector  technique.  To  weight  these  factors,  the 
following  equation  was  used: 

Ra(X)  - a • Rp(X)  ♦ 6 • RS(A)  ♦ y • RJA)  ♦ 6 • R$n(X) 

Where:  RA(A)  is  the  overall  local  spectral  albedo  reflectance, 

Rw(A)  is  the  average  spectral  reflectance  of  water, 

Rsn(A)  is  the  average  spectral  reflectance  of  snow,  and 
a.e.y.fi  are  fractions  of  each  surround  type  (a+B+y+fi  ■ 1.0). 

Figure  5-2  is  a typical  example  of  one  such  net,  weighted  spectral  reflectance 
of  the  terrain  surrounding  a measurement  site.  Spectral  reflectance  curves 
for  all  sites  can  be  found  in  Appendix  D.  Those  identified  as  "geo-botanical" 
were  reconstructed  using  red,  green,  and  blue  integrated  reflectances,  not 
from  the  photography,  but  from  maps  of  vegetation  and  ground  cover  for  the 
United  States. 

5.2  SURFACE  AEROSOL  OBSERVATIONS 
5.2.1  The  Albuquerque,  N.M.,  Experiment 

During  the  period  from  12  September  to  19  September  1976,  concurrent  ex- 
periments were  conducted  by  radiometric  laboratory  personnel  and  by  the  U.S. 
Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  N.M.  At 
the  request  of  the  A.F.  Wea  :her  Service,  the  White  Sands  Laboratory  provided 
particle  count  observations  to  both  parties  continuously  throughout  the 
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period  using  an  instrument  developed  by  Particle  Measurement  System,  Inc, 
called  a Classical  Scattering  Aerosol  Spectrometer.*  Charts  showing  data 
on  optical  scattering,  wind  speed,  and  relative  humidity  at  15  minute  inter- 
vals were  provided  to  this  contractor  through  AWS  direction.  These  three 
factors  tended  to  correlate,  but  showed  different  amounts  of  time  lag. 

Table  5-3  summarizes  these  observations  at,  or  just  preceding,  the  measure- 
ment period  for  the  four  days  that  both  experiments  were  conducted.  The 
average  of  the  morning  and  midday  RAOB  moisture  scale  heights  is  also  included 
in  Table  5-3. 

TABLE  5-3 

AEROSOL  OBSERVATIONS  AND  METEOROLOGICAL  DATA 


Date: 
Local  Time: 

9-12-76 

13:26-17:15 

9-13-76 

13:01-19:00 

9-14-76 

13-15-17:00 

9-19-7' 

13:17-19 

RAOB  Moisture  Scale  Ht.(cm) 

1.26 

1.30 

1.73 

1.39 

Avg  Particle  Concentration 
(Part/cm3) 

0.4 

1.30 

4.2 

1.5 

Peak  Particle  Concentration 
(Part/cm3) 

0.5 

4.3 

5.8 

2.3 

Avg.  Kind  Speed  (mph) 

9.0 

9.4 

7.4 

4.3 

Peak  Wind  Speed  (mph) 

15.0 

13.7 

9.0 

4.7 

Avg.  Wind  Speed  (mph)  during 
previous  24  hours 

— 

8.3 

14.1 

7.2 

Peak  Wind  Speed  (mph)  during 

— 

14.8 

26.6 

11.4 

previous  24  hours 


* Described  in  Research  and  Development  Technical  Report,  ECOM-5597, 
Response  Calculations  for  a Commercial  Light-Scattering  Aerosol 
Counter,  July  1976. 


5.2.2  Normalizing  Data  to  the  AFCRL  Aerosol  Model 


r 


To  normalize  the  particle  concentration  measurements  to  a definition  consistent 
with  the  AFCRL  model  atmosphere,  the  assumptions  of  the  particle  size  distribu- 
tion and  altitude  concentration  made  by  McClatchey  et  al.6  were  also  used  for 
our  observations.  Two  interpolations  were  required  before  particle  concentrations 
could  be  compared:  (1)  Although  the  particle  sizes  used  by  McClatchey  were  from 
.02  pm  to  10  pm,  the  measurement  sizes  used  at  Albuquerque  ranged  from  0.4  pm  to 
4 pm,  and  (2)  the  McClatchey  data  had  to  be  interpolated  for  a terrain  altitude 
of  1.5  km.  Although  factors  such  as  absolute  instrument  calibration  and  particle 
index-of-refraction  data  could  not  be  critically  compared,  these  two  first-order 
corrections  brought  the  measurements  and  theory  into  reasonable  agreement.  In 
Table  5-4,  data  from  McClatchey  can  be  compared  with  measured  data. 

TABLE  5-4 

COMPARISON  OF  NORMALIZED  PARTICLE  CONCENTRATION  DATA  AND  AFCRL  MODEL 


Date 

Measured  Mean 
Aerosol  Concentration 

Aerosol  Concentration 
Corrected  for  Size, 
Range,  and  Altitude 

Aerosol  Concentration 
AFCRL  Model 

(part /cubic  cm) 

(part/cubic  cm) 

(part/cubic  cm) 

9-12-76 

0.4 

317 

2828*  (Clear:  23  km 

9-13-76 

1.3 

1023 

visibility) 

9-14-76 

4.2 

3308 

13780**  (Hazy:  5 km 

9-19-76 

1.5 

1183 

visibility) 

As  can  be  seen  in  Tables  5-3  and  5-4,  there  is  a rough  correspondence  between 
particle  concentration  and  average  wind  speed  for  the  previous  24  hours  and  also 
possibly  with  moisture  scale  height.  Spectral  transmittance  data  for  these  same 


6 R.A.  McClatchey  et  al.  Optical  Properties  of  the  Atmosphere,  Third  Edition, 

Air  Force  Systems  Command,  AFCRL-72-0497,  1972. 

* Surface  concentration  equivalent  to  the  Clear  Model;  that  is,  23  km  visibility. 

Surface  concentration  equivalent  to  the  Hazy  Model;  that  is,  5 km  visibility. 
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four  days  showed  a marked  change  in  the  aerosol  coefficients  derived  from 
analysis  of  the  term  of  extinction  of  aerosol  scattering  for  these  days  (see 
section  2,  volume  3,  of  this  report).  A correction  for  the  aerosol  scattering 
extinction  coeficients  was  developed  on  the  basis  of  this  experiment. 


SECTION  6 


STORAGE  AND  RETRIEVAL  OF  COMPUTER  DATA 

6.1  DATA  TRANSFER  SOFTWARE 

It  was  not  possible  to  examine  or  analyze  spectral  and  map  data  collected 
by  the  mobile  radiometric  laboratory  in  the  cassette-tape  form  used  to 
transmit  data  from  the  field.  Therefore,  a series  of  programs  was  developed 
to  transfer  this  data  to  the  flexible  IBM  370  computer  system  and  to  manipu- 
late and  extract  this  data  once  it  had  been  stored  in  the  IBM  system. 

For  this  purpose,  data  in  terms  of  raw  measurement  amperes  and  in  engineer- 
ing units  were  transferred  from  the  cassette  tapes  to  nine-track  magnetic 
tapes  on  an  off-line  Data-General  computer  system  for  back-up  or  safe-storage 
purposes.  A second  transfer  was  also  made  of  the  data  to  a nine-track  mag- 
netic tape  in  engineering  units  that  could  easily  be  read  by  the  IBM  system. 
This  tape  was  then  used  to  prepare  an  IBM  3330  disk  pack  file  for  direct 
access  in  the  analysis  phases  of  the  project.  Figure  6-1  is  a schematic  of 
the  flow  of  data  after  the  data  was  received  at  this  contractor's  facility. 

6.2  DATA  GENERAL  COMPUTER  PROGRAM  STDG 

The  transfer  of  data  from  cassettes  to  nine-track  magnetic  tapes  for  back- 
up was  accomplished  by  the  program  STDG.  This  transfer  was  direct  with  no 
deletions  or  changes.  The  data  was  transferred  to  a 9-track  magnetic  tape 
using  an  off-line  Data  General  Computer  with  a 9-track  magnetic  tape  unit,  a 
Sykes  cassette  unit,  and  a teletype.  This  data,  in  amperes,  was  stored  on  a 
different  nine-track  tape  from  the  data  in  engineering  units.  Each  tape  was 
separated  into  files,  each  file  containing  one  day's  acquisitions.  These 
nine-track  tapes  were  an  exact  replicate  of  all  data  recorded  on  cassettes 
by  the  mobile  Radiometric  Laboratory.  Four  nine-track  tapes  were  used  to 
store  the  data:  two  tapes  for  data  in  amperes,  and  two  tapes  for  the  data 
in  engineering  units. 
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Figure  6-1.  Schematic  Showing  Flow  of  Data  in  Preparation  for  Analysis 
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6.3  DATA  GENERAL  COMPUTER  PROGRAMS  TAPDSK  AND  STIBM1 


To  do  large-scale  data  analysis,  it  was  necessary  to  transfer  the  data  in 
engineering  units  from  the  Data  General  computer  to  the  IBM  computer  (see 
flow  chart  in  Figure  6-1).  When  the  transfer  was  made,  an  extra  magnetic  tape 
was  made  in  addition  to  the  backup  tapes  for  two  reasons:  (1)  use  of  this 
extra  tape  permitted  the  backup  tapes  to  be  kept  in  storage  without  the  risk 
of  loss  of  its  data  from  manipulation  on  the  computer  system,  and  (2)  the 
format  of  the  backup  tapes  was  not  readily  readable  on  the  IBM  system  without 
a highly  specialized  assembler  program.  For  these  reasons,  it  was  more 
efficient  to  create  an  extra  tape  on  the  Data  General  computer  which  could 
read  these  tape  formats  directly. 

The  transfer  was  done  using  two  programs:  (1)  TAPDSK,  which  transferred  data 
from  the  nine-track  magnetic  tape  written  by  program  STDG  to  the  Data  General 
disk  file  STOOR,  and  (2)  STIBM1,  which  transferred  data  on  the  disk  file  STOOR 
to  nine-track  magnetic  tape  that  can  be  read  by  the  IBM  system.  The  data  was 
not  transferred  directly  from  cassette  to  nine-track  magnetic  tape  because  the 
Sykes  cassette  units  were  error  prone  and  often  required  special  handling  be- 
fore cassette  tapes  could  be  read  accurately. 

Whereas  the  program  TAPDSK  transferred  the  complete  file  from  the  nine-track 
magnetic  backup  tape  to  the  disk,  STIBM1  only  transferred  data  that  was  neces- 
sary for  analysis.  The  data  that  were  not  transferred  to  the  magnetic  tape 
were:  (1)  calibration  acquisitions,  (2)  acquisitions  that  were  prematurely 
terminated,  as  indicated  by  the  absence  of  a finish  time  in  the  header  block*, 
and  (3)  the  directory*  of  the  cassette. 
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STIBMl  changed  the  format  of  the  data  from  the  integer  format  to  floating 
point  format  (see  paragraph  4.5  of  Vol  1,  RELIB  Program).  The  program  also 
changed  the  comment  field  in  the  header  block  from  ASCII  to  EBCDIC.  These 
changes  made  it  possible  for  the  IBM  to  read  the  magnetic  tape  without  a new 

assembler  program. 

• 

6.4  COMPUTER  PROGRAM  STOVAN 

I 

■ 

6.4.1  Program  Description 

The  magnetic  tape  created  by  STIBMl  was  processed  on  an  IBM  370  System,  and 
the  data  was  stored  by  program  STOVAN  on  a disk  file  for  direct  access  in  the 
future.  STOVAN  stored  data  from  the  magnetic  tape  into  five  disk  files:  one 
file  for  each  of  the  four  instruments,  and  a file  for  the  sky  maps.  A sixth 
file  was  used  as  an  index  file  to  the  direct  access  files.  Other  functions 
performed  by  program  STOVAN  were  the  deletion  and  modification  of  data  already 
on  the  files,  the  computation  of  solar  altitude  and  solar  azimuth,  the  listing 
and  printer-plotting  of  data,  and  interpolation  of  the  data  from  10-nanometer 
to  5-nanometer  measurement  increments. 

6.4.2  Format  of  Data  Received  from  Program  STOVAN 

The  output  from  STOVAN  provides  files  in  three  formats:  (1)  the  index  file, 

(2)  the  sky  map  file,  and  (3)  the  fixed  geometry  files. 

The  format  for  the  Index  File*  consisted  of: 

Record  1 (unformatted)  : 

Array  52  by  25  (stored  as  integers  in  two-byte  words) : 

Contains  the  starting  record  number  for  each  week  of  the 
year  in  the  first  position  of  the  array  and  the  solar  alti- 
tude range  in  the  second  position  of  the  array  for  the  sky 
map  file. 

Next  (integer):  Record  number  of  the  next  unused  record 
in  the  sky  map  file. 


For  detailed  information  on  the  Index  File  Format,  see  IBM  System/ 
360  Operating  System  FORTRAN  IV  (G  and  H)  Programmer^  GuideV 
"Direct-Access-Programming." 


Record  2 (unformatted) : 


r 

I 
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Array  52  by  25  (integer  two  bytes) : Contains  a count  of  the 
number  of  records  for  each  week  of  the  year  in  the  first 
position  of  the  array  and  each  solar  altitude  range  in  the 
second  position  of  the  array. 

Records  3,  5,  7,  and  9 (unformatted): 

Each  record  corresponds  to  the  file  for  one  of  the  four 
instruments:  Sky  radiometer,  Irradiometer  1,  Irradiometer  2, 
and  Transmissometer. 

Array  12  by  20  by  25  (integer  two  bytes) : Contains  the  starting 
record  number  for  each  month  in  the  first  subscript,  a vacancy 
in  the  second  subscript  (not  used),  and  the  solar  altitude 
range  in  the  third  subscript. 

NEXT  (the  integer  variable):  Is  the  record  number  of  the  next 
unused  record  in  the  file. 

Records  4,  6,  8,  and  10  (unformatted) 

Each  record  corresponds  to  the  file  for  one  of  the  four  instru- 
ments: Sky  radiometer,  Irradiometer  1,  Irradiometer  2,  and 
Transmissometer . 

Array  12  by  20  by  25  (integer  two  bytes) : Contains  a count  of 
the  number  of  records  for  each  month  in  the  first  subscript  of 
array,  a vacancy  in  the  second  subscript  (not  used),  and  the 
solar  altitude  range  in  the  third  subscript. 





Each  record  on  the  sky  radiance  sap  data  file  has  the  following  format: 


Entry 

RVAN 

Byte 

(Integer:  2 bytes) 

Displacement  No. 

Description 

1-4 

Chain  (integer  4 bytes) 

Note 

1 

Record  No.  of  next 
record  in  series  of 
week  number  and  solar 
altitude  range. 

5,6 

Month 

5 

7,8 

Day 

6 

Date 

9,10 

Year 

7 

1 

11,12 

Hour 

8 

1 

GMT  at  beginning 

13,14 

Minute 

9 

of  acquisition 

15,16 

Second 

10 

17-20 

Latitude  (Real) 

Note 

2 

21-24 

Longitude  (Real) 

Note 

2 

In  degrees 

25-28 

Solar  Altitude  (Real) 

Note 

3 

29-32 

Solar  Azimuth  (Real) 

Note 

3 

33,34 

Wind  Speed 

38 

35,36 

Wind  Direction 

39 

37,38 

Humidity 

40 

Weather  Data 

39,40 

Temperature 

41  | 

41,42 

Pressure 

42 

43,44 

Elevation 

100 

Height  above  mean 
sea  level  (feet) 

45-204 

Comments 

102-181 

Up  to  160  characters. 

two  characters  per 
word. 

205,206 

First  Wavelength 

63 

207,208 

209,210 

Second  Wavelength 

Third  Wavelength 

64 

65 

Explicit  wavelengths, 
scan  mode  only 

Note  1:  For  detailed  information  on  the  Sky  Map  File,  see  IBM  System/360 
Operating  System,  FORTRAN  IV  (G  and  H)  Programmer  '"Guide,  •'Direct 
Access  Programming." 

Note  2:.  Program  STOVAN  changes  these  variables  from  the  RVAN  format  to  a more 
- usable  format. 

Note  3:  Program  RVAN  does  not  record  these  variables;  instead  they  are  com- 
puted internally  in  the  STOVAN  program. 
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Byte 

Entry 

(Integer:  2 bytes) 

RVAN 

Displacement  No. 

Description 

211-244 

(blanks) 

- 

- 

245-248 

Data  Point,  lrradiometer 

1 

249-252 

Data  Point,  lrradiometer 

2 

2S3-256 

Data  Point,  Transmissometer 

257-260 

261-264 

265-268 

Data  Point,  lrradiometer  1 

Data  Point,  lrradiometer  2 

Data  Point,  Transmissometer 

Note  4 

The  average  irradiance 
and  transmittance  mea- 
surements at  the  first, 
second,  and  third  wave- 

269-272 

Data  Point,  lrradiometer 

1 

lengths. 

273-276 

Data  Point,  lrradiometer 

2 

277-280 

Data  Point,  Transmissometer 

> 

281-474 

(blanks) 

- 

- 

475-476 

Azimuth 

20 

First  Instrument 

477-478 

Azimuth 

23 

Second  Instrument 

479-480 

Azimuth 

26 

Third  Instrument 

481,482 

Azimuth 

29 

Fourth  Instrumment 

483,484 

(blanks) 

- 

- 

485,486 

Elevation 

21 

First  Instrument 

487,488 

Elevation 

24 

Second  Instrument 

489,490 

Elevation 

27 

Third  Instrument 

491,492 

Elevation 

30 

Fourth  Instrument 

493,494 

(blanks) 

- 

- 

495,496 

No.  of  Points 

50 

Number  of  points  in 
full  map  table;  scan 
mode  only. 

497,500 

(blanks) 

- 

- 

501,2036* 

Map  Table 

Note  4 

Array  of  sky  radiance 
data  (After  all  data 
for  one  wavelength  is 
recorded,  recording  of 
the  next  wavelength  is 
done  without  a break.) 

2037-2054 

(blanks) 

- 

- 

Note  4:  These  variables  have  no  corresponding  RVAN  header  block  displacement 
numbers,  but  they  can  be  found  in  data  files. 

* See  Vol  1,  Para  4.3  for  the  order  of  the  data. 
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Each  record  in  the  spectral  data  files  of  the  Sky  Radiometer,  Irradiometer  1, 
Irradiometer  2,  and  Transmissometer  has  the  following  format  for  fixed  geometry 
acquisitions: 


Entry  RVAN 


Byte 

(Integer:  2 bytes) 

Displacement 

No. 

Description 

1-4 

Chain 

Note 

1 

The  record  No.  of  the 
next  record  in  the 
series  organized  by 
month  and  solar  alti- 
tude range. 

5,6 

Month 

5 

7.8 

Day 

6 

Date 

9,10 

Year 

7 

11,12 

Hour 

8 

13,14 

Minute 

9 

GMT  at  start  of  acqui- 

sition 

15,16 

Second 

10 

17-20 

Solar  Altitude 

Note 

3 

21-24 

Blank 

25-28 

Latitude 

Note 

2 

In  degrees 

29-32 

Longitude 

Note 

2 

In  degrees 

33-36 

Solar  Azimuth 

Note 

3 

37-38 

Instrument  Azimuth 

20, 23, 26, or 

29 

39,40 

Instrument  Elevation 

21, 24, 27, or 

30 

Note  1:  For  detailed  information  on  the  Sky  Map  File,  see  IBM  System/360 
Operating  System,  FORTRAN  IV  (G  and  H)  Programmer’s  Guide,  "Direct 
Access  Programming." 

Note  2:  Program  STOVAN  changes  these  variables  from  the  RVAN  format  to  a 
more  usable  format. 

Note  3:  Program  RVAN  does  not  record  these  variables ; instead  they  are 
computed  internally  in  the  STOVAN  program. 
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Entry 

(Integer:  2 bytes) 

RVAN 

Displacement  No.  Description 

41,42 

Nind  Speed 

38  , 

43,44 

Wind  Direction 

39 

| Weather  Data 

45,46 

Humidity 

40 

47,48 

Temperature 

41  ' 

| 

49,50 

Pressure 

42 

f 

51,52 

Sea  Elevation 

100 

Height  above  mean 
sea  level  (feet) 

53-212 

Comments 

102-181 

Up  to  160  characters 
2 characters  per 
word 

213-216 

(blanks) 

- 

- 

217-864 

Radiance,  Irradiance,  or 
Transmittance 

Note  4 

Data  in  5 nm  incre- 
ments. 

865-876 

(blanks) 

- 

- 

6.5  COMPUTER  PROGRAM  EXTRACT 
6.5.1  Program  Description 

After  the  data  is  stored  in  the  IBM  370  Computer,  the  program  EXTRACT  is  used 
to  retrieve  it.  The  main  purpose  of  the  program  is  to  extract  data  by  month, 
day,  year,  solar  altitude,  time,  and  instrument  type,  thus  providing  data  for 
only  that  specific  measurement  set.  Alternatively,  a complete  listing  of  data 
for  all  dates  and  times  can  be  obtained  by  specifying  only  the  irradiance  type 
(ex:  "DH",  for  daylight  horizontal).  The  program  also  lists  and  printer-plots 
data  on  command.  The  EXTRACT  program  reads  the  disk  files  created  by  STOVAN, 
extracts  the  requested  data,  and  writes  it  in  a second  sequential  file  in  the 
same  format  as  STOVAN  but  without  the  variable  termed  CHAIN.  Use  of  the  second 
file  permits  more  rapid  access  to  data- than  would  be  possible  by  providing 
access  to  the  original  STOVAN  disk  file.  This  was  the  normal  procedure  for  the 
analysis  phase  of  this  project  where  speed  in  handling  the  large  volume  of  data 
was  of  importance. 


Note  4.  These  variables  have  no  corresponding  RVAN  header  block  displace 
ment  numbers,  but  they  can  be  found  in  the  RVAN  data  files. 
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6.5.2  EXTRACT  Program 

Because  this  program  has  not  been  documented  elsewhere,  its  input  require- 
ments have  been  included  in  this  report  for  purposes  of  completeness. 

6. 5. 2.1  CARD  No.  1 Input.  Col.  1-4:  PLOT,  or  blanks.  If  "PLOT"  is  speci- 
fied, all  the  data  that  is  extracted  is  also  listed  and  plotted;  otherwise 
the  data  is  extracted  only  and  written  to  the  new  disk  file. 

6. 5. 2. 2.  CARD  No.  2 Input.  Only  columns  73-75  must  be  specified.  If  nothing 
is  entered  in  columns  1-72,  all  data  corresponding  to  the  specified  instrument 
will  be  extracted.  If  the  data  corresponding  to  a range  of  dates  is  desired, 
columns  1-6  are  the  starting  date  and  columns  7-12  are  the  ending  date.  If 
no  range  is  desired,  columns  7-12  are  left  blank.  Similarly,  if  data  is 
desired  for  a range  of  times,  columns  13-18  are  the  beginning  hour,  minute, 
and  second;  and  columns  19-24  are  the  ending  hour,  minute,  and  second.  If 
no  range  is  desired,  columns  19-24  are  left  blank.  Data  for  a range  of  solar 
altitudes  can  also  be  extracted  if  columns  25-31  specify  the  starting  solar 
altitude  and  columns  33-39  specify  the  ending  solar  altitude.  If  no  range  is 
desired,  columns  33-39  are  left  blank.  If  specific  ranges  are  desired,  the 
following  format  for  data  inputs  is  followed: 


Col. 

1-2: 

Starting  month  (right  justified) 

Col. 

3-4: 

Starting  day  (right  justified) 

Col. 

5-6: 

Starting  year 

Col. 

7-8: 

Ending  month,  or  blank  (right  justified) 

Col. 

9-10: 

Ending  day,  or  blank  (right  justified) 

Col. 

11-12: 

Ending  year 

Col. 

13-14: 

Starting  hour  (right  justified) 

Col. 

15-16: 

Starting  minute  (right  justified) 

Col. 

17-18: 

Starting  second  (right  justified) 

Col. 

19-20: 

Ending  hour,  or  blank  (right  justified) 

Col. 

21-22: 

Ending  minute,  or  blank  (right  justified) 

Col. 

23-24: 

Ending  second,  or  blank  (right  justified) 

Col. 

25-31: 

Solar  altitude 

Col. 

32:  (blank) 
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Col.  33-39:  Solar  altitude,  or  blank 

Col.  40-56:  (blanks) 

Col.  57-63:  Latitude  (If  data  for  a particular  site  is  required, 

the  latitude  should  be  specified.) 

Col.  64:  (blank) 

Col.  65-72:  Longitude  (If  data  for  a particular  site  is  required, 
specify  longitude.) 

Col.  73-75:  SKY  (Sky  radiance  fixed  geometry  mode) 

IR1  (Irradiometer  1 data) 

IR2  (Irradiometer  2 data) 

IRA  (Transmissometer  data) 

MAP  (Sky  radiance  map  mode) 

Col.  76:  (blank) 

Col.  77-79:  DH  (Daylight  horizontal  for  IR1  or  IR2) 

DV  (Daylight  vertical  for  IR1  or  IR2) 

SH  (Skylight  horizontal  for  IR1  or  IR2) 

SVB  (Skylight  vertical  backlit  for  IR1  or  IR2) 

SVF  (Skylight  vertical  frontlit  for  IR1  or  IR2) 


APPENDIX  A 
PARAMETER  SUMMARY 


The  following  table  gives  the  exact  values  or  ranges  by  day-site 
for  the  parameters  important  to  this  project. 


TABLE  A-l 


PARAMETER  SUMMARY 


SITE 

DAT 

EB 

■riM 

TERRAIN 

xefle:tav:e 

• 65c yy. 

(TTriV/Hi 

Little  lock,  Ark. 

8-6-75 

2005  to  2245 

58.73  to  27.62 

1.989 

2.874 

3.440 

311 

0.085 

Athens,  Ga. 

8-12-75 

1803  to  2336 

70.33  to  8.76 

2.978 

- 

3.457 

811 

0.085 

Athens,  Ca. 

8-13-75 

2201  to  2215 

28.24  to  25.22 

2.965 

3.784 

3.256 

811 

0.080 

Athens,  Ca. 

8-14-75 

1342  to  1608 

33-56  to  61.80 

2.698 

3.008 

3.786 

811 

0.080 

Tampa,  Fla. 

8-19-75 

1429  to  2215 

74.11  to  23-06 

4.169 

4.576 

3.808 

20 

0.078 

Apalachicola,  Fla. 

8-27-75 

1447  to  2400 

70.34  to  1.13 

3.883 

4.319 

3.447 

20 

0.078 

Apalachicola,  Fla. 

8-28-75 

2015  to  2340 

48.92  to  5.02 

4.009 

4.112 

3-874 

20 

0.078 

Denver,  Col. 

10-21-75 

1855  to  2355 

39.45  to  2.31 

0.302 

0.405 

0.593 

5400 

0.123 

Denver,  Col. 

10-22-75 

1915  to  0003 

38.66  to  O.83 

0.338 

- 

0.940 

5400 

0.123 

Tint  Hill,  Ta. 

11-6-75 

1751  to  2005 

33-57  to  19.37 

1.823 

2.452 

2.522 

420 

0.080 

Vint  Hill,  Ta. 

11-9-75 

1615  to  2145 

34.25  to  2.30 

1.759 

2.145 

2.173 

420 

0.080 

Tint  Hill,  Ta. 

11-11-75 

1520  to  2120 

33.71  to  6.29 

0.626 

0.6’1 

1.073 

420 

0.080 

St.  Cloud,  Kin. 

2-21-76 

1715  to  2345 

33-70  to  0.17 

0.312 

0.338 

O.386 

1028 

0.407 

St.  Cloud,  Min. 

2-22-76 

1952  to  2345 

33-70  to  Cu43 

0.427 

0.261 

0.405 

1028 

0.407 

St.  Cloud,  Kin. 

2-23-76 

1415  to  2000 

34.43  to  10.33 

0.749 

0.713 

0.718 

1028 

0.407 

St.  Cloud,  Kin. 

2-24-76 

1734  to  2345 

34.80  to  0.80 

0.762 

0.433 

0.510 

1028 

0.407 

St.  Cloud,  Kin. 

2-27-76 

1650  to  1845 

39.39  to  28.94 

I.O64 

- 

0.557 

1028 

0.407 

Glasgow,  Mont. 

3-4-76 

1415  to  2000 

4-59  to  35.57 

0.242 

0.202 

0.283 

2279 

0.407 

Glasgow,  Mont. 

>5-76 

1400  to  1909 

2.44  to  35.92 

0.221 

0.230 

0.578 

2279 

0.407 

CuiU ajTJte,  Hash. 

4->76 

1706  to  2015 

30.62  to  47.60 

1.091 

0.994 

0.877 

181 

0.038 

Quillsyute,  Hash. 

4-5-76 

1532  to  1805 

16.50  to  39.29 

0.860 

0.981 

1.3a 

181 

0.038 

Quillayute,  Hash. 

4-6-76 

2300  to  0145 

10.93  to  36.81 

1.101 

0.852 

0.772 

181 

0.038 

Alameda,  Calif. 

4-l>76 

1533  to  2240 

22.30  to  61.53 

1.071 

1.095 

0.920 

19 

0.112 

Alameda,  Calif. 

4-20-76 

1349  to  1430 

3.57  to  11.38 

0.931 

2X33 

0.968 

19 

0.112 

Alameda,  Calif. 

4-21-76 

1350  to  1645  ' 

. 3.78  to  38.04 

0.793 

0.752 

1.251 

19 

0.112 

Alameda,  Calif. 

4-26-76 

1326  to  0255 

0.22  to  65.96 

0.530 

0.472 

0.567 

19 

0.112 

Vandenberg,  Calif. 

>7-76 

1930  to  0025 

28.57  to  72.27 

1.957 

1.442 

1.639 

328 

0.324 

Vandenberg,  Calif. 

>12-76 

1343  to  0255 

O.64  to  73.57 

1.600 

1.804 

1.171 

328 

0.324 

San  Diego,  Calif. 

>19-76 

2050  to  0020 

70.53  to  28.02 

1.221 

1.386 

1.153 

408 

0.145 

San  Diego,  Calif. 

>26-76 

1935  to  0025 

27.80  to  78.44 

1.021 

1.460 

1.871 

408 

0.145 

San  Diego,  Calif. 

6-1-76 

1832  to  2015 

70.50  to  77.62 

1.133 

1.481 

1.5U 

408 

0.145 

Tucson,  Axis. 

6-7-76 

2041  to  0225 

70.48  to  0.18 

1.400 

1.320 

1.080 

2705 

0.263 

Tucson,  Axis. 

6-8-76 

1515  to  0130 

10.49  to  80.73 

0.988 

0.768 

1.290 

2705 

0.263 

Tucson,  Axis. 

6->76 

1536  to  2Q30 

39.34  to  80.81 

1.022 

1-352 

1.079 

2705 

0.263 

Tucson,  Aria. 

6-11-76 

1226  to  1930 

0.95  to  80.96 

0.871 

0.645 

0.652 

2705 

0.263 

Tucson,  Aria. 

6-17-76 

1221  to  1905 

0.00  to  80.24 

1-385 

1.507 

1.924 

2705 

0.263 

Vinslow,  Aria. 

6-28-76 

1215  to  1930 

0.10  to  78.24 

1.296 

1.400 

1-612 

4895 

0.232 

TABLE  A- 1 (CONT'D) 


Little  Rock,  Ark. 

7-14-76 

1412  to  1627 

35.8  to  63.OO 

B9 

Little  Rock,  Ark. 

7-19-76 

1120  to  1620 

1.27  to  61.12 

ESS 

Little  Rock,  Ark. 

7-20-76 

1128  to  1240 

3.00  to  18.60 

3.426 

Little  Rock,  Ark. 

7-21-76 

1112  to  1800 

-0.4  to  75.07 

2.918 

Athene,  Ga. 

8-4-76 

1110  to  1457 

3.74  to  50.06 

1.999 

Athens,  Oa. 

8-6—76 

1111  to  1515 

3.61  to  53.30 

2.643 

Athens,  Ga. 

8-10-76 

1430  to  1705 

44.38  to  70.16 

2.849 

Athens,  Oa. 

•JJ-76 

1100  to  1600 

0.73  to  60.90 

2.737 

Athens,  Ga. 

8-12-76 

1055  to  1350 

-O.43  to  35.36 

2.777 

Athens,  Ga. 

8-13-76 

1055  to  1600 

-0.60  to  60.46 

3-433 

Apalachicola,  Fla. 

8-22-76 

1225  to  2355 

3.12  to  59.00 

3-981 

Apalachicola,  Fla. 

8-23-76 

1245  to  1725 

18. 90  to  70.99 

3.891 

Albuo.uerque,  N.K. 

9-12-76 

1326  to  1715 

6.67  to  49.81 

1.403 

Albuquerque,  N.V.. 

9-13-76 

1301  to  1900 

1.55  to  58.44 

1.097 

Albuquerque,  N.K. 

9-14-76 

1315  to  1700 

4-44  to  47.32 

1.737 

Albuquerque,  K.K. 

9-19-76 

1317  to  1938 

4.15  to  56.13 

1.385 

Ogden,  Uteh 

9-27-76 

1607  to  1920 

28.84  to  46.91 

1.252 

Ogden,  Utah 

9-28-76 

1336  to  0120 

-1.83  to  46.52 

1.013 

Ogden,  Utah 

9-2^76 

1330  to  2330 

0.35  to  46.13 

0.894 

Dayton,  Ohio 

10-22-76 

1200  to  1530 

0.06  to  32.70 

0.421 

Dayton,  Ohio 

10-26-76 

1422  to  1515 

23-19  to  29.94 

Dayton,  Ohio 

10-28-76 

1205  to  2230 

-0.19  to  36.75 

Dayton,  Ohio 

10-29-76 

1430  to  1630 

23.44  to  35.18 

0.789 

r^y  • 
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APPENDIX  B 


SURFACE  METEOROLOGICAL  DATA 


The  following  table  summarizes  the  surface  meteorological  observations 
nearest  the  time  of  the  Rawinsonde  launches  of  1200,  1800,  and  2400 
GMT.  Values  that  are  asterisked  are  for  observations  taken  within 
hours  before  or  after  the  time  of  launch;  however,  fifty  six  percent  of 
these  observations  were  within  1 hour  of  launch.  All  other  values  were 
essentially  concurrent  with  the  RAOB  launches. 
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TABLE  B-l 

SURFACE  METEOROLOGICAL  DATA 


tat* 

Location 

Surfaca 

Rraaauia 

(■b) 

Ta^aratura 

(5) 

12002  15002  00002 

Ralativa  Hiaidity 
(5) 

12002  18002  00002 

-TT — 

Mind 

Spaad 

(Knott) 

Avg.  Kind 
Diraction 
(Dagr.  from  T.N.) 

5/12/75 

Athana 

954 

91 

89* 

48 

52* 

1.45 

60 

1/14/75 

Ath*ns 

961 

54* 

90* 

61* 

49* 

4.0 

20/240 

«/ 19/75 

Taapa 

954 

94 

95* 

48 

51* 

12.47 

500 

5/27/75 

Apalachicola 

9S0 

91* 

91 

90* 

56* 

61 

55* 

7.17 

120 

5/25/75 

Apalachicola 

980 

89* 

57* 

5.0 

100 

10/21/75 

D*nv*r 

769 

82* 

80* 

21* 

19* 

1.14 

— 

10/22/75 

D*nv*r 

755 

82 

83* 

21* 

20* 

1.76 

... 

11/9/75 

Vint  Hill 

955 

79 

54 

1.56 

160 

11/11/75 

Vint  Hill 

9S1 

65 

27 

1.17 

200 

2/21/76 

St.  Cloud 

903 

25 

28* 

61 

S6* 

6.S 

340 

2/22/76 

St.  Cloud 

906 

50* 

55* 

57* 

50* 

0.45 

20 

2/25/76 

St.  Cloud 

591 

52* 

45 

71* 

44 

2.47 

240 

2/24/76 

St.  Cloud 

884 

52 

49* 

SO 

SS* 

0.0 

... 

5/4/76 

Glasgou 

844 

-6* 

4* 

71* 

61* 

0.85 

320 

5/5/76 

Glasgow 

541 

2* 

12* 

75* 

66* 

1.42 

140 

4/S/76 

Quillayuta 

866 

SS 

SS 

1.95 

180 

4/6/76 

Quillayuta 

897 

60 

38 

1.96 

20/220 

4/21/76 

A 1 seeds 

1014 

57* 

66* 

2. 25 

20 

4/26/76 

A i seeds 

inis 

60* 

64 

72 

36* 

31 

25 

10.12 

20 

5/7/76 

Vsndenberg 

1003 

65* 

63 

57* 

67 

7.00 

300 

S/12/76 

Vandanbarg 

1003 

71* 

8S 

82 

51* 

56 

40 

1.96 

20 

5/19/76 

San  Oitgo 

998 

66 

62 

7.04 

20 

5/26/76 

San  Diego 

1003 

80* 

42* 

3.60 

20 

6/7/76 

Tucson 

920 

101* 

102 

18* 

17 

6.56 

20 

6/5/76 

Tucson 

921 

91* 

94 

100 

22* 

19 

18 

0.78 

20 

6/9/76 

Tucson 

923 

92 

24 

2.41 

20 

6/11/76 

Tucson 

922 

66* 

80 

52* 

25 

0.50 

20 

6/17/76 

Tucson 

920 

51* 

95 

21* 

22 

5.74 

20 

6/25/76 

Vinslow 

854 

72* 

96 

25* 

20 

0.56 

20 

7/14/76 

Littl*  Rock 

997 

88* 

94* 

61* 

55* 

0.17 

260 

7/19/76 

Littla  Rock 

1009 

84* 

52* 

0.79 

60/140 

7/20/76 

Littl a Rock 

1009 

78 

77 

0.0 

— 

7/21/76 

Littla  Rock 

1005 

79 

76 

90 

2.76 

260 

5/4/76 

Athant 

994 

69 

78 

2.51 

40 

5/6/76 

Athens 

990 

72 

74 

1.15 

220 

5/10/76 

Athens 

991 

79* 

58* 

59* 

51* 

0.75 

20 

5/11/76 

Athens 

995 

7S 

85* 

76 

60* 

2.54 

60 

5/12/76 

Athant 

99S 

75 

81 

0.12 

100 

5/15/76 

Athens 

989 

76 

89* 

78 

51* 

0.26 

240 

5/22/76 

Apalachicola 

1016 

90* 

91 

49* 

54 

5.55 

260 

5/25/76 

Apalachicola 

1015 

81 

96 

71 

34 

1.98 

340 

9/12/76 

Albuquerque 

559 

71* 

74* 

40* 

34# 

6.97 

20 

9/15/76 

Albuquerque 

840 

62* 

85 

46* 

25 

1.90 

20 

9/14/76 

Albuquerque 

544 

75* 

82* 

55* 

55* 

0.0 

— 

9/19/76 

Albuquerque 

555 

70* 

51 

44* 

25 

1.10 

240 

9/27/76 

Ogdan 

571 

69 

59 

10.49 

100 

9/25/76 

Ogden 

572 

55* 

70 

75 

S2* 

55 

25 

5.99 

100/520 

9/29/76 

Ogdan 

575 

54* 

71 

78* 

48* 

51 

25* 

5.05 

100/500 

10/26/76 

Dayton 

990 

44* 

70* 

6.10 

20 

10/28?76 

Dayton 

990 

30 

44 

47* 

77 

58 

54* 

2.55 

220 

10/29/76 

Dayton 

992 

45* 

45* 

49* 

44* 

5.55 

220 

Raadings  ut  within  2.5  hour!  of  ala  tiao 


APPENDIX  C 


MEASUREMENT  INDEX 


This  index  lists  all  the  spectroradiometric  measurements  made  by  the 
mobile  radiometric  laboratory  during  the  collection  phase  of  the  project. 
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OATE  TIME  INST  s*  OU*L  COMMENTS 
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APPENDIX  D 

SURROUND  SPECTRAL  REFLECTANCES 


Graphs  D-l  through  D-14  show  spectral  reflectances  of  the  terrain  surround- 
ing the  measurement  site,  as  estimated  from  U-2  color  photography  or  geo- 
botanical  maps. 


Figure  D-3.  Local  Terrain  Spectral  Albedo  Reflectance  for  Appalachicola,  Fla. 
and  Ta*pa,  Florida  (Geo-Botanical) 
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Terrain  Spectral  Albedo  Reflectance  for  Denver,  Colorado 
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Figure  D-6.  Local  Terrain  Spectral  Albedo  Reflectance  for  Quillayute,  Washington 


Figure  D-9.  Local  Terrain  Spectral  Albedo  Reflectance  for  San  Diego,  California 
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